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Extensions of Einstein gravity with higher-order derivative terms arise in string theory and other 
effective theories, as well as being of interest in their own right. In this paper we study static 
black-hole solutions in the example of Einstein gravity with additional quadratic curvature terms. 
A Lichnerowicz-type theorem simplifies the analysis by establishing that they must have vanishing 
Ricci scalar curvature. By numerical methods we then demonstrate the existence of further black- 
hole solutions over and above the Schwarzschild solution. We discuss some of their thermodynamic 
properties, and show that they obey the first law of thermodynamics. 


The well-known problem of the non-renormalisability 
of Einstein gravity has given rise to many attempts to 
view it as an effective low-energy theory that will re¬ 
ceive higher-order corrections that become important as 
the energy scale increases (see, for example, m)- In 
string theory, the Einstein-Hilbert action is just the first 
term in an infinite series of gravitational corrections built 
from powers of the curvature tensor and its derivatives. 
In other approaches, only a finite number of additional 
terms might be added. It was shown in [5] that if one 
adds all possible quadratic curvature invariants to the 
usual Einstein-Hilbert action one obtains a renormalis- 
able theory, albeit at the price of introducing ghost-like 
modes. Arguments have been given for why these might 
not be fatal to the theory (for example, see [5] for a recent 
discussion). In any case, it is worthwhile to study in de¬ 
tail the properties of the theory of Einstein gravity with 
added quadratic curvature terms, in order shed light on 
the question of whether it has irredeemable pathologies 
or whether they can be controlled in some manner. 

Black holes are the most fundamental objects in a 
theory of gravity, and they provide powerful probes 
for studying some of the more subtle global aspects 
of the theory. It is therefore of considerable interest 
to investigate the structure of black-hole solutions in 
theories of gravity with higher-order curvature terms. 
In this paper, we report on some investigations of 
the static, spherically-symmetric black-hole solutions 
in four-dimensional Einstein-Hilbert gravity with added 
quadratic curvature terms, for which the most general 
action can be taken to be 

1 = J (yi? - , (1) 

where a, /3 and 7 are constants and is the Weyl 


tensor. We shall work in units where we set 7 = 1, and 
the equations of motion following from Q are then 

Rfii/ -\- ‘IjSRi,Rpi/ i^RUn^) 

+2/3(gp,ni? - VpV,i?) = 0 , (2) 

where B^i, = (V^V'^ -b \R^'^)Cppu (7 is the Bach tensor, 
which is tracefree. 

In general, the theory describes a system with a mas¬ 
sive spin-2 mode with mass-squared m| = 1/ ( 20 ?) and a 
massive spin-0 mode with mass-squared ttiq = l/(6/3), in 
addition to the massless spin-2 graviton. These massive 
modes will be associated with rising and falling Yukawa 
type behaviour in the metric modes near infinity [?], of 
the form and In particular, one can 

expect that if generic initial data is set at some small 
distance, the rising exponentials will eventually domi¬ 
nate, leading to singular asymptotic behaviour. In seek¬ 
ing black-hole solutions, the question then arises as to 
whether the rising exponentials can be avoided for ap¬ 
propriately finely-tuned initial data. 

It can easily be seen that any solution of pure Einstein 
gravity will also be a solution of ([^ , and so in particular 
the usual Schwarzschild black hole continues to be a so¬ 
lution in the higher-order theory. The question we wish 
to address, then, is whether there exist any other static 
black hole solutions, over and above the Schwarzschild 
solution. 

Static, spherically-symmetric black-hole solutions have 
been investigated in , using generalisations of the Lich- 
nerowicz and Israel theorems for Einstein gravity. Since 
we will arrive at somewhat different conclusions, we shall 
briefly summarise the key elements in [5], although de¬ 
rived in a different notation. We consider static met¬ 
rics of the form ds\ = —\^dt^ -\- hijdx^dx^, where A 
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and hij are functions only of the three spatial coordi¬ 
nates a;h Taking the trace of the field equations ([^ gives 
/3 (□ — mo)-R = 0. We then multiply this by \R and inte¬ 
grate over the spatial domain from a putative horizon out 
to infinity. Expressed in terms of the covariant derivative 
Di with respect to the spatial 3-metric hij, this gives 


In our work, we have therefore carried out a numerical 
investigation of the solutions. To do this, we begin by 
supposing that there exists a black-hole horizon at some 
radius r = Tq > 0, at which the metric functions h and 
/ vanish, and we then obtain near-horizon Taylor expan¬ 
sions for h{r) and f{r), of the form 


Vhd^a 


D\\RD,R) - \{D,Rf - ml\R^ =0. (3) 


Since A vanishes on the horizon, it follows that if DiR 
goes to zero sufficiently rapidly at spatial infinity the to¬ 
tal derivative {i.e. surface term) gives no contribution, 
and the non-positivity of the remaining terms then im¬ 
plies i? = 0. In other words, as shown in 0 , any 
static black-hole solution of Q must have vanishing Ricci 
scalar. This leads to a great simplification, and it means 
that one can, without loss of generality, study the case 
of pure Einstein-Weyl gravity {i.e. Q with /3 = 0), since 
obviously the term quadratic in R makes no contribution 
to the held equations for a conhguration with R = 0. Fur¬ 
thermore, the trace of the held equations ([^ for Einstein- 
Weyl gravity immediately implies i? = 0. In fact, the two 
diherential equations for h and / are both now of only 
second order in derivatives. 

The second stage of the discussion in 0 then involved 
looking at the remaining content of ([^, i.e. the non¬ 
trace part. According to 0 , this led to another integral 
identity that then implied, under certain assumptions, 
that R^i, = 0. If this were correct, then the conclusion 
would be that the usual Schwarzschild solution is the 
only static black hole solution of the theory described by 
0- However, we hnd that there are sign errors in the 
expression given in 0 . Setting i? = 0, as already argued 
above, multiplying ([^ by XR^'', and then integrating 
over the spatial region outside the horizon gives 


Vhd^ 


- \X{D,R - 4D^Rijf + 4X{D^R,jf 


-4X{D[,Rj]kf + X{D,Rjkf - \XR^{ml + R) 
-X{mlKiR,j - ^R^^jkR^i)] = 0 , 


( 4 ) 


where Wi = XR^^DiRjk + \XRDiR — 2XR^^DjRik — 
XRDjRi^, and R is the Ricci scalar of the spatial metric 
hij. Although the surface term will give zero, the mix of 
positive and negative signs in the bulk terms precludes 
one from obtaining any kind of vanishing theorem for 
the Ricci tensor of the four-dimensional metric. This 
raises the intriguing possibility that there might in fact 
exist static, spherically symmetric black-hole solutions 
over and above the Schwarzschild solution. 

The equations of motion following from Q are too 
complicated to be able to solve explicitly, even for the 
case of the static, spherically-symmetric ansatz 

dr 2 

= —h{r)dt^ + -Tj-T + r‘^{d9‘^ + sin^ Odcj)'^). (5) 

f{r) 


h{r) 
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{r - To) -k ha (r - -k hg (r - rof 


/(^) = fi{r- ro) + f 2 {r- + fs {r - rof H-(6) 


Substituting into the equations of motion ([^, with /3 set 
to zero for the reasons discussed above, the coefficients 
hi and fi for i > 2 can be solved for in terms of the 
two non-trivial free parameters rg and fi. There is also 
a “trivial” parameter, corresponding to the freedom to 
rescale the time coordinate, which we have accordingly 
written in the form of an overall scaling of h{r). Thus we 
have 


1 - 2/i rp I - /i rp _ I-2/irg 3(1-/irp) 
firl Sa/frp ’ ^ rg 8a/i rp 


and so on. (We used Taylor expansions to 0{{r — rp)®) 
in our numerical integrations.) The Schwarzschild solu¬ 
tion corresponds to fi = 1/rp, and so it is convenient to 
parameterise fi as 


fi 


1 + d 

ro 


( 7 ) 


with non-vanishing S characterising the extent to which 
the near-horizon solution deviates from Schwarzschild. 

We use the expansions to set initial data at a radius r^ 
just outside the horizon, and then use numerical routines 
in Mathematica to integrate the equations out to large 
radius. Generically, one finds that for a given choice of 
the parameters rp and S the solution rapidly becomes 
singular as one integrates outwards from r = ri, as ex¬ 
pected in view of our earlier observations about the rising 
Yukawa terms in the asymptotic form for the metric. If 
we fix a particular value for rp, we can then use the shoot¬ 
ing method to try to home in on a special value of S for 
which the outward integration can proceed without en¬ 
countering a singularity. Of course in practice, because 
of accuracy limitations in the integrations, the solution 
will always eventually become singular at large enough 
r. The signal for a good black-hole solution is that /(r) 
and h{r) should approach constants as r increases, (in 
fact Ricci scalar flatness implies /(r) must approach 1), 
and that by stepping up the accuracy and precision goals 
in the calculations one can extend at will the maximum 
upper limit r = rf for which the smooth behaviour can 
be achieved. In practice, by running the routines with 
accuracy and precision goals of order 20 decimal places, 
we have been able to obtain very clean and trustworthy 
solutions out to at least 60 times the horizon radius. 

Our findings are that there exists a range of values for 
the horizon radius, bounded below by a certain multiple 
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of the length \fa^ for which we can obtain precisely one 
static black-hole solution in addition to the Schwarzschild 
solution. In order to make the statement of our results in 
the most concise possible way, it is convenient, without 
loss of mathematical generality, to make a specific choice 
for the value of a in Q. We shall take 

^=\- 

We then find that for each choice of rg > r™'", where 

« 0.876 , (9) 


M 




FIG. 1: The masses as functions of vq (left plot), and 
as functions of the Hawking temperatures (right plot) for 
the Schwarzschild (dashed line) and non-Schwarzschild (solid 
line) black holes. 


we can find a non-Schwarzschild static black hole. For 
each such rg, there is a corresponding value d = d* of 
the “non-Schwarzschild parameter” that yields the black- 
hole solution without a singularity at spatial infinity. As 
rg is taken closer and closer to the value r™’", the re¬ 
quired value S* becomes smaller and smaller, tending 
to zero at rg = Thus the Schwarzschild and the 

non-Schwarzschild black holes “coalesce” as rg = r™™ is 
approached. 

As rg is increased above r™™, the Schwarzschild and 
non-Schwarzschild black-hole solutions separate more 
from one another (and in particular the required value 
of S increases). The mass of the Schwarzschild black hole 
is simply ^rg, and thus it increases linearly as rg in¬ 
creases. By constrast, the mass of the non-Schwarzschild 
black hole decreases as rg increases, until at rg = r™=° 
it becomes massless, where 

rS'=°« 1.143. (10) 

(The definition of mass in higher-derivative theories was 
discussed in li HJ. For asymptotically-flat black holes 
it is just i the coefficient of 1 /r in gtt (assuming gu 
is normalised canonically at infinity).) Interestingly, 
if rg is increased beyond r“^°, one can still obtain a 
non-Schwarzschild black-hole solution for an appropri¬ 
ate choice of S, but now the mass is actually negative. 
In other words there is still a regular horizon, and the 
metric is asymptotically flat at large distances, but the 
metric function / now rises above 1 as r increases from 
Ti, before sinking down to 1 again in the asymptotic re¬ 
gion. Figure 1 shows the masses of the Schwarzschild 
and non-Schwarzschild black holes as a function of rg, 
and their masses as functions of their Hawking tem¬ 
peratures. The maximum possible mass for the non- 
Schwarzschild black hole, attained when rg = r™“, is 
given by = ^r™'" « 0.438. From the slope of 

M{T) it can be seen that the specific heat C = dM/dT 
is negative for both black holes, and more negative for the 
non-Schwarzschild black hole at a given temperature. 

The plots of the metric functions / and h for the ex¬ 
amples of a positive-mass black hole with rg = 1 , and 
a negative-mass black hole, with rg = 2 , are shown in 
Figure 2. 



FIG. 2: The non-Schwarzschild black hole for rg = 1 (left 
plot) and rg = 2 (right plot). In each plot the upper curve is 
f(r) and the lower curve is h{r). For clarity we have chosen a 
rescaling of h so that it approaches |, rather than 1, to avoid 
an asymptotic overlap of the curves. 


Having established the existence of the non- 
Schwarzschild black holes, it is instructive to study some 
of their thermodynamic properties, and to compare these 
with the properties of the Schwarzschild black holes. In 
order to do this, we have collected the numerical re¬ 
sults for a sequence of black-hole solutions with rg in 
the range r“™ « 0.876 < rg < 1.5, and then fitted 
the data to appropriate polynomials. Because we are 
working with a higher-derivative theory, the entropy is 
not simply given by one quarter of the area of the event 
horizon, and instead we need to use the formula derived 
by Wald O |9]. This has been evaluated for the ansatz 
(§ in quadratic curvature gravities in [ 101 , applied 
to our case with /? = 0 and 7 = 1 in Q this gives 
S = 7rrg -I- 47ra(l — /i rg) = 7rrg — 47rai5*. (There is a 
freedom to add a constant multiple of the Gauss-Bonnet 
invariant to the Lagrangian, which shifts the entropy by 
a parameter-independent constant without affecting the 
equations of motion. We have used this to ensure the 
entropy of the Schwarzschild black hole vanishes when 
the mass vanishes.) We then find that the mass and the 
temperature of these non-Schwarzschild black holes, as a 
function of the entropy, take the form 

M « 0.168-f 0.131 S'- 0.00749 S'2- 0.000139 , 

T « 0.131 - 0.0151 S - 0.000428 S^ -f • • • . (11) 

It can be seen that dM/dS » 0.131 — 0.0150 S — 
0.000417 S^, which is very close to the expression for the 
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temperature. Thus the non-Schwarzschild black holes are 
seen to obey the first law cLM = TdS to quite a high preci¬ 
sion. Note that the expressions for M and T as a function 
of S for the Schwarzschild black holes are very different 
in form, with M = ( 5 '/ 47 r)i /2 and T = i( 7 rS')-i/ 2 . 


s r' 



FIG. 3: The first plot shows the entropy as a function of 
mass, and the second shows the free energy F = M — TS 
as a function ofT, for the Schwarzschild (dashed line) and 
non-Schwarzschild (solid line) black holes. 

It is interesting to note that the entropy of the non- 
Schwarzschild black hole of a given mass is always less 
than the entropy of the Schwarzschild black hole of the 
same mass. The two entropies approach each other 
asymptotically as rg approaches r™'". This can be seen 
in the left-hand plot in Figure 3. It is also of interest 
to look at the free energy F = M — TS as a function 
of temperature. This is shown in the right-hand plot in 
Figure 3. It can be seen that the free energy is always 
larger for the non-Schwarzschild black hole at a given 
temperature, with the two curves again meeting at the 
lower limit when rg = 

In this paper, we have used black holes to probe some 
of the consequences of interpreting the action 0 as com¬ 
plete classical action in its own right. We have seen that 
there exists a second branch of static, spherically sym¬ 
metric black holes, over and above the Schwarzschild so¬ 
lutions. These are not Ricci flat, although they do have 
vanishing Ricci scalar. Restoring the factors of a and 7 
that we hxed in our numerical simulations, the second 
branch of black holes have masses, which can become 
negative, bounded approximately by M < 0 . 438 -^/ 207 . 
Thus in a regime where a is small, which one might 
hope would correspond to a small correction to Einstein 
gravity, the second branch of black holes will be tiny, 
and will actually have very large curvature near the hori¬ 
zon, thus tending to invalidate the requirement that the 
curvature-squared should be small. The fact that their 
mass can be negative, violating the usual positive-mass 
theorem of standard Einstein gravity, indicates that the 
ghost-like nature of the quadratically-corrected action is 
becoming dominant in this regime; one may view the 
negative-mass black holes as condensates dominated by 
the contribution of the ghost-like massive spin -2 modes. 
It could be viewed as a satisfactory outcome of our in¬ 
vestigation that the only indications of the existence of 
black holes with potentially pathological properties in the 


quadratic-curvature theories occur in a regime where yet 
higher-order corrections, as in string theory, are going 
to be important also. It would be interesting to obtain 
analytical proofs of the existence of the numerical solu¬ 
tions we have found. Although this could be challenging, 
it might perhaps be easier to obtain restricted no-hair 
theorems that confirmed the apparent absence of non- 
Schwarzschild black holes outside the parameter range 
where we have found them. 

Naturally, the fourth-order equations of motion ([^ 
have wider classes of solutions than the black-hole solu¬ 
tions with horizons that we have considered here. These 
were initially investigated in [ 1 ] and will be given a more 
detailed analysis, along with a more extensive analysis of 
the black hole solutions, in m 
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